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Abstract 


rhis document contains graphical data necessai) lot flic pie!timnai\ design of ballistic 
missions to Jupitei. ('ontouis of launch energv icquiremenfs. as well as mam othei 
launch and .lupitci amval paiameteis. ate piesented m launch date ai rival date space 
toi a!l launch oppoilunities fioni 1 1 >N> tluough 2005. In addition, an extensive text is 
included which explains mission design methods, tiom launch window development to 
Jupitei piolv and oihitei amval design, utilizing the giaplucal data m this volume as 
well as numciiHis equations lelating vaiioiis paiameteis. This is the fust of a planned 
seiies M mission design documents which will appl> to all planets and some othei bodies 
in tlte solai s\stem. 



Preface 


This publication is one of a series of volumes devoted to interplanetary trajectories of 
different types. Volume I, of which the present publication is Part 3, describes ballistic 
trajectories. Parts 1. 2, and 4, which will be published in 1983. will treat ballistic trajec- 
tories to Venus. Mars, and Saturn, respectively. 
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I. Introduction 

The purpose of tins series of Mission Design Handbooks is 
to provide trajectory designers and mission planners witli 
graphical trajectory information* sufficient for preliminary 
interplanetary mission design and evaluation. In most respects 
the series is a continuation of the previous three volumes of 
the Mission Design Data. TM 33-736 (Ref. 1 ) and its predeces- 
sors (e.g., Ref. 3); it extends their coverage to departures 
through the year 3.003 A.D. 

The entire series is planned as a sequence of volumes, each 
describing a distinct mission mode as follows: 

Volume 1- Ballistic (i.e.. unpowered) transfers be- 
tween Karth and a planet, consisting of 
one leg trajectory arcs. For Venus and 
Mars missions the planet-tu-Furth return 
trajectory data are also provided. 

Volume II: Gravity-Assist (G'A) trajectory transfers. 

comprising from two to four ballistic 
interplanetary legs, connected by suc- 
cessive planetary swinghys. 

Volume III: Delta- V-EGA (AVEGA) transfer trajec- 
tories utilizing an impulsive deep-space 
phasing and shaping burn, followed b> a 
return to Earth for a G/A swingbv maneu- 
ver taking the spacecraft (S O to the 
eventual target planet. 

Each volume consists of several parts, descnbing trajecto: 
opportunities for missions toward specific taigei ot swingby 
bodies. 

This Volume 1. Pait 3 of the series is devoted to ballistic 
transfers between Earth and Jupiter. It ucscnhcs tiajectories 
taking from 1 to 5 yeais of flight time for the 30 successive 
mission opportunities, departing Earth in the following years: 
|‘>K5, 1%. IWl 1 4 >KK. 1W, | W(). IWj. IU‘P'3. 1W3 4. 
1004 5. |Wi>, JW7. iws.iwo :t)0(\ 3001. 3003. 3003 4. 
3004 5. and 3005 u 

Individual vaiiahles piesented herein are described in detail 
in subsequent sections and suininari/ed again in Section IV. 
Suffice it to say here that all the data aie presented in sets of 
11 contour plots each, displayed on the launch date aiuval 
date space foi each opportunity, kequired departure cncigy 
C x% departure asvinpiot. decimal ion and light ascension, 
aiuval 1^ and its equatonal directions, as well Sun and 
Eaith direction angles with iespect to the depart me ainxal 
asvmpioies. are presented. 

It should he •* ied that pan the lauiuh %pace cor, ..d 
may require launcher energies not presently available. 


but certainly not unrealistic using future orbital assembly 
fwhnuuu 1 *. 

a 

A separate series of volumes (Ret, 3) is being published 
concurrently to provide purely geometrical (i.e., trajectory- 
independent) data oi planetary positions and vi?wing/oricnta- 
tion angles, experienced by a spacecraft in the vicinity of these 
planetary bodies. The data cover the time span through 3030 
A.D.. in order to allow sufficient mission duration time lor all 
Earth departures, up to 3005 A.D. 

The geometric data are presented in graphical form and 
consist of 30 quantities, combined .. * > eight plots for each 
calendar year and each target planet. Tiu /raphs display equa- 
torial declination and right ascension of Earth and Sun (plan- 
etocentric). as w ell as those of the target planet (geocentric); 
heliocentric (ecliptic) longitude of the planet, its heliocentric 
and geocentric distance; cone angles of Earth and Canopus, 
clock angle of Earth (when Sun,C anopus-oric: i ted); Earth-Sun- 
planel. as well as Sun-Earth-planet angles; and finally, rise and 
set times for six deep-space tracking stations assuming a h-deg 
horizon mask. This information is similar to that in the second 
part of each of the volumes previously published (Ref. 1 ). 

II. Computational Algorithms 

A. General Description 

The plots for the entire series were computer-generated. A 
minimum of editorial and graphic suppoit was postulated from 
the outset in an effort to reduce cost. 

A number of computer programs were created and/or modi- 
fied to suit the needs of the Handbook production. 

The computing effort involved the generation ot arrays ot 
transfer trajectory arcs connecting departure and arrival 
plane is on a large number ot suitable dates at each body. 
Algorithms t computational models) to solve this problem can 
vary greatly as to their complexity, cost ot data generated, and 
resulting data accuracy. In light of these considerations, the 
choice of methods used m this etfort has been assessed. 

B. Two-Body Conic Transfer 

h ich departure anoal date combination sepreseiiis a unique 
transfer tiaiectoiy between two qv ,itied bodies, it the mi.ubei 
ot (evolutions ot the space*, tat t about the primal v te e . the 
Sun) is s|H\ Med I he I ainbeit theorem pi ovules a suitable 
tiamewoik foi the imputation ot siulr tiaict tones, hut 
only it restricted two bodv wonu motion presaiU 

Restricted two-hodv motion implies that tire dynamical 
system consists ot onlv two bodies, one ot which, tire pnuiaiv. 


1 
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is so much more massive than the other, that :i!! of the system’s 
gravitational attraction may he .ssumed as concentrated at a 
point the center of t li .at primary body, the secondary body 
of negligible mass (e.g., the spacecraft) then moves in Kep- 
lerian (conic) orbits about the primary (e.g.. the Sun) in such 
a way that the center of the primary is located at one of the 
foci of the conic (an ellipse, paiabola, or hyperbola). 

The Lambert Theorem states that given a value of the ; *uvi- 
tatinnai parameter id (also known as Ghf) for the central body, 
the time of tlight between two arbitrary points in space, 
Rj and R 2 , is a function of only three independent variables: 
the sum of the distances of the two points from the focus, 
|R , 1 + IR 2 I, the distance between the two points C= |R 2 - R, 1. 
and the semimajor axis, a . of the conic orbital flight path 
between them (Fig. 1.) 

Detailed algorithm descriptions of the Lambert method, in* 
eluding necessau branching and singularity precautions, are 
presented in numerous publications, e.g.. Refs. 2 and 4. The 
computations result in a set of conic classical elements [a, t\ /. 
il, gj. v x ) and the transfer angle. Au l2 . or two equivalent 
spacecraft heliocentric velocity vectors. V;, v one at depar- 
ture. the other at the ainval planet. Suhtiaction of th«* ippio- 
priate planetary heliocentric velocity vector. V* . at 

the two corresponding times from each of these two space- 
craft velocity vectors results in a pair of planetocentric velocity 
states "at infinity'" with respect to each planet (Fig. 2) 

V ~. I ' W,, Hi.A*n . ( ” 

where i = 1 and 2 refer to positions at departure and arrival, 
respectively. The scalar of this L^ vector is also referred to as 



Fig. 1. The Lambert problem geometry 


the hyperbolic excess velocity, ” F*in Unity" or simply "speed” 
(e.g.. Ref. 4). The represents the velocity of the spacecraft 
at a great distance from the planet (where its gravitationa 1 
attraction is practically negligible). It is attained when the 
spacecraft has climbed away from the departure planet, fol- 
lowing injection at velocity \ t : 



or before it starts its fall into the arrival planet’s gravity well, 
where it eventually reaches a closest approach < periapse, C/A ) 
velocity F /f . 


y 

S IK Ilf 



Fig. 2 . Departure geometry and velocity vector diagram 
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1 ho variables r t and r p refer to the departure injection and 
arrival periapse plane tocen trie radii, respectively. Values for 
the gravitational parameter /i (or GM) are given in subsequent 
Section V on constants. 

The V*. vectors, computed by the Lambert method, repre- 
sent a body center to body center transicr. They can, however, 
be translated parallel to themselves at either body without 
excessive error due to the offset, and a great variety cf rt .istic 
departure and arrival trajectories may thus be constructed 
through their use. to be discussed later, fhe magnitude and 
direction of (’„ as we!l as the angles that this vector forms 
with the Sun and Larth direction vectors at each terminus, are 
required for these mission design exercises. 

Miso * the relatively small terrestrial planets are suited 
to 1. ' analyzed by the Lambcr t method, as the problem can be 
adequately represented by the restricted two-body formula- 
tion, resulting in flight time errors of less than I day an 
accuracy that cannot even b*r read from the contour plots 
presented in this document. 


C. Pseudostate Method 

Actual precision interplanetai transfer traiectories. espe- 
cially those involving the giant outer planets, do noticeably 
violate the assumptions inherent in the Lambert Theorem. 
The restricted two-body problem, on which that theorem is 
based, is supposed to describe the conic motion of a massless 
M.vonJar\ (i.e.. the spacecraft) about the point mass of a pri- 
mary attractive body (i.e.. the Sun), both objects being placed 
in an otherwise empt) Universe. In reality, the gravitational 
a* 'ractmn ot either departure or tatget body may significantly 
ah-.r the entire transfer trajectory. 

Numerical V-hody truicctory integration could be called 
upon to represent the true physical mode! tor the laws of 
motion, but would be too costly, lonsidci lug the number ot 
complete Iijjcs. tones requited to fully search and describe a 
given mission opportunity 

The pseudostate theory, lust introduced by S. Wilson 
(Kef. 5) and modified to solve the tluee-body Lambert prob- 
lem by l). V By ines (Ret e>). represents an extremely useful 
improvement over the standard (.amber t solution. Lot the 
giant planet missions, it can correc about percent tit the 
three-body errors incurred, e g., up to 30 days in flight tune on 
a Jupitci -bound tourney 


Pseudostate theory is based on the assumption that for 
modest gravitational perturbations the spacecraft conic motion 
about the primary and the pscudo-conic displacement due to a 
third body may be superimposed, if certain rules are followed. 

The method, as applied to transfer trajectory generation, 
does not provide a flight path only its end states. It solves the 
original Lambert problem, however, not between the true 
planetary positions themselves, but instead, between two com- 
puted "pseudostates.*' These are obtained by iteration on two 
displacement vectors off the planetary ephemeris positions on 
the dates *»f departure and arrival. By a suitable superposition 
with a planetocenfnc rectilinear impact hypetbola and a 
constant-velocity. "zero gravity," sweep back at each end of 
the Lambertian conic (see Fig. 3). a satisfactory match is 
obtained. 


Of the five arcs involved in the iteration, the last three 
(towards and at Jupiter) act over the lull flight tune. 2^7 J2 . 
and represent: 


(1) Conic heliocentric motion between the two pseudo- 
states R ( and R* (capital R is used here for all helit 
centric positions). 

(2) The transformation of R* to a planctocentric position, 
r 2 (lower case r is used for planetocentric positions), 
performed in the usual manner is followed by a 
"constant-velocity" sweepback in time to a point 
r* * r 2 - X &T l2 . correcting the pi anetocen ‘ im- 
position r 2 to what it would have been at Tj. had 
there been no solar attraction during ^7^. and 
finally . 


(3) The planetocentric rectilinear incoming hyperbola, 
characterized by incoming ('-infinity V^, .t radial 



r* RECTILINEAR 


HYPERBOLA. 


/ JUPITER 


Z r 2 ^2 ^2 


Fig. 3. Psaudostata transfer gaomatry 
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target planet impact, and a trip time A7' |2 from r J 
to perijove. which can be satisfied by iteration on the 
r 2 -magnitude and ;hus also on R*. 


This last aspect provides for a great sim plit1c.it ion of the 
formulation as the R? end-point locus now moves only along 
the || rl vector direction. The resulting reduction in com- 
puting cost *s significant, and the equivalence to the Lambert 
point-to-point conic transfer model is attractive. 


The first two segments of the transfer associated with the 
departure planet may be treated in a like manner. If the planet 
is Larth. the pseudostate correction may he disregarded fie . 
Rj = R | ). or else the duration of Earth's perturbative effect 
may be reduced to a fraction of It can also be set to 

equal a fixed quantity, e g., ^7j = 20 days. The latter value 
was in fact used at the Earth's side ot the liansfe* in the 
data generation process for this document 


The rectilinear pseudostatc method thus involves an itera 
five procedure, utilizing the standard Lambert algorithm to 
obtain a starting set of values for at each end of the trans- 
fer arc This first guess is then improved by allowing the 
planetocentric pseudostate position vector r* to be scaled up 
and down, using a suitable partial at either bodv. such that 
A T r . the :imc required t*> fall along the rectilinear hypci- 
bola through rj (the sum of the sweepback distance l\ mt 
X A7j and the planetocentric distance ir* |). equal the giavita- 
tional perturbation duration. AT Both r* and V along 
which the rectilinear fall occuis. are continuously icset utiliz* 
” i the latest values of magnitude and direction of at each 
end. i. ot the new Lambert transfc aic. as the iteration pn> 
gresses. The proceduie converges rap'dly as the hyperbolic trip 
time discrepancy. AA/' = A - A T . tails below a preset 
small tolerance 


Once the V. vectors at each planet are converged upon, 
the desired output variables can be generated and contour 
plotted by existing standard algorithms. 


III. Treiectory Characteristics 
A. Mission Space 

All realistic launch and inaction vehicles are energy limned 
and impose very stringent constraint on the interplan.ta’v 
mission selection process. Only those transfer opp. 
which occur near the times of a minimum Larth . •* 
energy requ.icmeni are thus of practical interest. O* ■ -•Si*-* 
side of such an optimal date, departure energy mere.-. *s. first 
slowly, followed bv a rapid increase, thus requiring eitiier a 
greatei launch capability, or alternatively a lower allowable 
pavload mass. A "launch period.** measured in days or even 


weeks, is thus definable: on any day within its confines the 
capability of a given launch/injection vehicle must equal or 
exceed the departure energy requirement for a specified 
payload weight. 


In the course of time these minimum departure energy 
opportunities do recur regularly, at “synodic period ’ intervals, 
reflecting a repetition of the relative angular geometry of the 
two planets. If u>, and lj 2 arc the orbital angular rates of the 
inner and outer of the two planets, respectively, moving about 
the Sun in circular orbits, then the mutual configuiation of 
the two bodies changes at the following rate: 

Lj i2 = w | " . rad/s 

If a period of revolution. P , is defined as 


then 

J_ _ J_ JL 

P s !\ /\ 

where P s . the synodic period, is the period of planetary 
geometry recurrence, while P x and P 2 are the orbital “sidereal 
(i.e., inertial) periods'* of the inner (faster) and the outer 
(slower) planet considered, respectively. 

Since planetary orbits are neither exactly circular nor co- 
planar, launch opportur : 'ies do not repeat exactly, some 
years heing better than others in energy requirements or in 
other parameters. A complete repeat of trajectory character- 
istics occurs only when exactly the same orbital geometry of 
departure and arrival body recurs. For negligibiy perturbed 
planets approximately identical inertial positions in space at 
departure and arrival imply near-recurrence of transfer tra- 
jectorv characteristics. Such events c^n rigorously be assessed 
onlv for nearly resonant nonprcccssing planet;* r y orbits, i.e.. 
tor those whose periods can he related in terms of integer 
fiactions. For instance, if five revolutions of one body cor- 
respond to three revolutions of the i ther, that time interval 
would constitute the “period of repeated characteristics. 14 
Near-integer ratios provide nearly petitive configurations 
with respect to the lines ot apsides and nodes. The Earth- 
relative synodic period ot Jupiter is b>S 4 days, i.e., about 
U months. Lach cycle ot 11 consecutive Jovian mission 
opportunities amounts to 44 k 7 7 d.*vs and is nearly repetitive, 
driven bv Jupiter’s sidereal period ot 1 I StOO veals u 4442 b 
days It is obvmus that tor an identical mutual angulai georn 
etiv Jupiter would he tuund beyond its uurtial posit r in m the 
previous vVi.ie bv l days worth ot motion (4 vs dee I wtu!** 


(4) 


(51 


(h) 
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Faith would have com pie led 4387.7/365.25 - 12.013 icvolu 
Hons, being ahead ol ihe old maik by the same angulat aiiiounl 
asiupitei 

A variety ol* consideiatioiis lore 1 the realist t launch period 
not to occui at the ininimnni energy comhm >tion ol depai 
turc and arrival dates, launch vehicle readiness status, proce- 
dure slippage, wearier anomalies, multiple launch strategics, 
arrival characteristics all cause the launch or. more generally, 
the departure period to be extended ovei a number of days or 
weeks and not necessarily centered on the minimum energy 
date. 

For this document, a 1 20-day departure date coverage span 
was selected, primarily in order to encompass launch energy 
requirements ol* up to a = 200 km*/* 2 contour, where 
C\ = rj i.e., twice the injection energy per unit mass. 


/: | l i 12. Arrival date coverage was set at 1600 days to dis- 
play missions t'lom 1- to 5-years flight time. 

TIip matrix of departure and arrival dates to he presented 
comprises the “mission space” lor each departure opportunity. 

B. Transfer Trajectory 

As previously stated, each pair of departure/arrival dates 
specifies a unique transfer trajectory. Fach such point in the 
mission space has associated with it an array of descriptive 
variables. Departure cnctgy. charactcri/i J by is by far the 
most significant among these parameters. It increases towards 
the '*dges of the mission space, hut it also experiences a dra- 
matic rise along a “ridge “ passing diagonally from lower left 
to upper right across the mission space (Fig 4). This dislur- 


1AMTH JUPITiH UflMy .in 
HA! 1 1 ST 1 C TRAIMSrtH THAJICHUtV 
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fig. 4. Mission space in deparlure/ar rival date coordinates, typical example 
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bance is associated with all diametric, i.e., near-IKO-deg, out of the ecliptic and forces a polar 180-deg transfer, in order 
transfer trajectories (Fig. 5). to pick up the target’s vertical out-of-plane displacement. 


In 3-uimensional space the fact that all planetary orbits are 
not strictly coplanar causes such diametric transfer arcs to 
require high ecliptic inclinations, culminating in a polar flight 
path for an exact 180-deg ecliptic longitude increment be- 
tween departure and arrival points. The reason for this 
behavior is, as shown in Fig. 5, that the Sun and both trajec- 
tory end points must lie in a single plane, while they are also 
lining up along the same diameter across the ecliptic. The 
slightest target planet orbital inclination causes a deviation 

= TRUE TRANSFER 
ANGLE 

fa = TRANSFER ANGLE 
POLAR IN ECLIPTIC PLANE 



Fig. 5. Effect of transfer angle upon inclination of 
trajectory arc 
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Fig. 6. Nodal transfer geometry 


The obvious sole exception to this rule is the nodal transfer 
mission, where depaiture occurs at one node of the target 
planet orbit plane with the ecliptic, whereas arrival occurs at 
the opposite such node. In these special cases, which recur 
every half of the repeatability cycle, discussed in the preceding 
paragraph, the transfer trajectory plane is indeterminate and 
may as well lie in the departure planet’s orbit plane, thus 
requiring a hsser departure energy (Fig. 6). The opposite strat- 
egy (i.e., a transfer in the arrival planet's orbit plane) may be 
preferred if arrival energy, V„ 2 . is to be minimized (Fig. 7). 

It should be noted that nodal transfers, being associates 
with a specific Jupiter arrival date, may show up in the data on 
several consecutive Jovian mission opportunity graphs (e.g., 
1991-1994/5). at points corresponding to the particular nodal 
arrival date. Their mission space position moves, from oppor- 
tunity to opportunity, along the 180-deg transfer ridge, by 
gradually sliding towards shorter trip times and earlier relative 
departure dates. Only one of these opportunities would occur 
at or near the minimum departure energy or the minimum 
arrival V w date, which require a near-perihelion to near- 
aphelion transfer trajectory. These pseudo-Hohmann nodal 
transfer opportunities provide significant energy advantages, 
but represent singularities, i.e.. single-time-point missions, 
with extremely high error sensitivities. Present-day mission 
planning does not allow single fixed-time departure strategies; 
however, future operations modes, e.g.. space station “on- 
time" launch, or alternately Earth gravity assist (repeated) 
encounter at a specific time, may allow the advantages of a 
nodal transfer to be utilized in full. 

The 180-deg transfer ridge subdivides the mission space 
into two basic regions: the Type I trajectory space below the 
ridge, exhibiting less th?n 180-deg transfer arcs, and the 
Type II space whose transfers are longer than 180 deg. In 
general the first type also provides shorter trip times. 

Trajectories of both types are further subdivided in two 
parts-Classes 1 and 2. These are separated, generally hori- 
zontally, by a boundary representing the locus of lowest 
C 3 energy for each departure date. Classes separate longer 
duration missions from shorter ones within each type. Type I, 
Class 1 missions could thus be preferred because of their 
shorter trip times. 

Transfer energies become extremely high for very short 
trip times, infinite if launch date equals arrival date, and of 
course, meaningless for negative trip times. 

The reason that high-inclination transfers, as found along 
the ridge, also require such high energy expenditures at depar- 
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Fig. 7. Mission space with nodal transfer 


ture is that the spacecraft velocity vector due to the Earth's 
orbital velocity must be rotated through large angles out of the 
ecliptic in addition to the need to acquire the required transfer 
trajectory energy. The value of C 3 on the ridge is large but 
finite; its saddle point minimum value occurs for a pseudo- 
Hohmann (i.c., perihelion to aphelion) polar transfer, requiring 

C 3 = I'l +lj * 2400 km 2 /s 2 (7) 

where 1^ = 2 < >.7ftb km/s, the Earth's heliocentric orbital 
velocity, and a n = 5.2 AU. Jupiter's (the arrival planet's) semi- 
majoi axis B> a siniilat estimate, it can be shown that for a 
true nodal pseudo-Mohniann tiansfei. the iniminum eneigy 
uv * d would reduce to 



This is the lowest value of C 3 required to fly from Earth to 
Jupiter, assuming circular planetary orbits. 


Arrival t-infinity. V«,^, is at its lowest when the transfer 
trajectory is near-coplanar and tangential to the target planet 
orbit at arrival. 


Both C 3 and near the ridge can be significantly low- 
ered if deep-space deterministic maneuvers are introduced into 
the mission. Tne “broken-plane" maneuvers are a category of 
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such ridge-counteracting measures, which can nearly eliminate 
all vestiges ot the near-180-deg transfer difficulties. 

The basic principle employed in broken-plane transfers is to 
avoid high ecliptic inclinations of the trajectory by performing 
a nlane change maneuver in the general vicinity of the halfway 
point, such that it would correct the spacecraft's aim toward 
the target planet’s out-of-ecliptic position (Fig. 8). 

Graphical data can be presented for this type of mission, 
but it requires an optimization of the sum of critical AV 
expenditures. The decision on which AVs should be included 
must be based on some knowledge of overall staging and arrival 
intentions, e.g., departure injection and arrival orbit insertion 
vehicle capabilities and geometric constraints or objectives 
contemplated. As an illustration, a sketch of resulting contours 
of C 3 is shown in Fig. 9 for a typical broken-plane oppor- 
tunity represented as a narrow strip covering the ridge area 
(Class 2 of Type I and Class 1 of Type II) on a nominal 1985 
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Fig. 8. Broken-plane transfer geometry 
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launch/arrival date C 3 contour plot. The deep-space maneuver 
&V BP is transformed into a C 3 equivalent by converting the 
new broken-plane C 3flp to an injection velocity at parking 
orbit altitude, adding and converting the sum to a new 

and slightly larger C 3 value at each point on the strip. 

C. Launch/Injection Geometry 

The primary problem in departure trajectory design is to 
match the mission-required outgoing K-infinity vector, V^, to 
the specified launch site location on the rotating Earth. The 
site is defined by its geocentric latitude, 0 7 . and geographic 
east longitude, X L (Figs. 10 and 1 1). 

Range safety considerations prohibit overflight of popu- 
lated or coastal areas by the ascending launch vehicle. For each 
launch site (e.g.. Kennedy Space Center, Western Test Range, 
or Guiana Space Center), a sector of allowed azimuth firing 
directions is defined (measured in the site's local horizontal 
plane, clockwise from north). For each launch vehicle, the 
allowed sector may be further constrained by other safety 
considerations, such as spent stage impact locations down the 
range and/or down-range significant event tracking capabilities. 

The outgoing K-infinity vector is a slowly varying function 
of departure and arrival date and may be considered constant 
for a given day of launch. It is usually specified by its energy 


magnitude C 3 = |VJ 2 , called out as C 3 L in the plots and 
representing twice the kinetic energy (per kilogram of injected 
mass) which must be matched by launch vehicle capabilities, 
and the K-infinity direction with respect to the inertial Earth 
Mean Equator and equinox of 1950.0 (EME50) coordinate 
system: the declination (i.e., latitude) of the outgoing asymp- 
tote 6^ (called DLA), and its right ascension (i.e., equatorial 
east longitude from vernal equinox, T) (or RLA). These 
three quantities are contour-plotted in the handbook data pre- 
sented in this volume. 

1. Launch azimuth problem. The first requirement to be 
met by the trajectory analyst is to establish the orientation of 
the ascent trajectory plane (Ref. 7). In its simplest form this 
plane must contain the outgoing V-infinity (DLA, RLA) 
vector , the center of Earth , and the launch site at lift-off 
(Fig. 10). As the launch site partakes in the sidereal rotation of 
the Earth, the continuously changing ascent plane manifests 
itself in a monotonic increase of the launch azimuth. with 
lift-off time, /} , (or its angular counterpart,^ - a L , measured 
in the equator plane): 

cos 0 7 X tan 5^ - sin X cos (o^ - a L ) 


(9) 



Fig. 10. Launch/injection trajectory plane geometry 
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Hie dailx tunc lustoix of a/unuth can he obtained t '10111 Fq. 
(*>) for a given a^, departure asymptote direction by cate- 
lu IK following quadrant rules explained below and using the 
follow ing au.xiliaiy expressions (see Fig 1 1 ) 




(»7M /m /7 =■- 100 075542(10 + 0>>S5M7 J4h0 X 


i 2 . 1 UU 5 \ 10 ,J \ i/:„ 
>o 


(ID 


wheie 


u f - light ascension of launch site (0 ; , \ ; ) at / ; 
(deg) 

(i 7 /.*I /m n . - Greenwich hour angle at l/ 1 lIMTof the 
launch date, the eastward angle between 
vernal equinox and the (iteenwich meridian 
(deg) 

*/ MI = launch date in terms of integer da\s since 
O' 1 Jan. I. !°>0 (da\ s) 

I ant suleieal rotation late of Faith ( 15.05 lOoo 1 * 
deg h tor tire l'*Ss-2lH)'' lime period) 

htt oil tune (h. lAI 1 ) 


A relative launch time, / K// . measured with lespect to an 
inertial refer ence (the depaituie asxmptote mendiariX right 
ascension u_ ), can be defined as 


K1 t 


24 0 - 


u', 


. Il 


(12) 


/ I Kill 


This tune represents a generalized sidereal time of launch 
elapsed since the launch vile last passed the depaituie asymp 
tote meridian, u fc 


Ihe actual liuvnwich Mean lime (GMl) of launch, : t , 
max be oht amed tiom t Kl r bx adding a date. site, and asxmp 
tote dependent adiustment 


\i M.vn 

! ! *K/i * 

/ IK/// 

Ihe expression lor (lq 'M must be used with computa 
tional regard tin quadrants. singular points, and sign comen* 
Mom It the launch is known to be direct tie. eastwaull. 
then 1, , when cot an is uegatixe. must be collected to 

1. , ♦ IS00 lor letiogiade (westward) launches. 

/ \ j i- 

IS0 0 dee must be added in the thud (when cotan 2 is post 
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LAUNCH SITE 
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Fig. 11. Earth equator plane definition ot angles involved 
in the launch problem 


live) and otiO.O deg in the fourth quadrant (when cotan 1 
is negative). 

A generalized plot ot relative launch time f K/ r vs bunch 
azimuth 2^ can be constructed based on Fqs. t 1 *) and (12). 
if a fixed launch site latitude is adopted (e.g., 2N.5 deg 

lot Kenned x Space Might 1'entci at Cape Canaveral. Florida), 
Such a plot is picscutcd in Fig. 12 with depaituie asxmptote 
declination as the contour paiatuetei. The plot is applicable 
to any leahstic departure condition, independent of update, 
or Hue launch tune /, (Ret S). 

2. Daily launch window:. Inspection of Fig. 12 indicates 
that genet all} two contours exist tot each declination value 
(e.g.. -10 deg), one occunmg at r during the a in. 

hours, the othei in the p m. bouts of the asxmptote lelatne 
’Max." 

Since lilt-oil tunes are bounded bx preselected launch site 
dependent limiting values ot launch a/unuth 2.', (e g , "0 deg 
and 1 15 deg), each of the two declination contouis thus con- 
tains a segment duitng which launch is pemussible "a launch 
window.** Flic two segments on the plot do define the two 
ax ail able daily launch window %. 

\s can be seen liom ft ig 12. loi X >% . 0. the two dailx 

launch oppoitumties ate separated bx exactlx 12 hours, with 
an increasing ! the\ close in on each otliei. until at i tS ^ 1 
,0/ thex merge into a single daib oppoiluml} ft ; oi |»S J • 
\0 t . a ■‘spirt'* ot that single launch xxmdow occuis. disallow 
mg an ever increasing sector ot a/unuth values Hus sector is 
sxmmettic about east and Us limits can be detci mined tiom 
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Fig. 12. Generalized relative launch time y vs launch a/iniuth and departure asymptote declination iS^ 

Pair of typical example launch windows tor iS^, 10 deg are shown by bold cuive segment 

(reproduced from Ref. 8). 

v ci ' s ^ «• The angular equatorial distance between the ascending no** ’ 

sM1 / i\tn cos O * a,u ^ *hi- launch site meridian is given h\ 


As gets longer, the sector of unavailable launch azimuths 
reaches the safely boundaries ol permissible launches, and 
planar launch ceases to exist (Fig. l.M. This subject will be 
addressed again in the discussion of “dogleg** ascents. 

Figure 14 is a sketch ol a typical daily launch geometry 
situation, shown upon a Mercator map ol the celestial sphere. 
The two launch windows exhibit a similar geometiv since the 
inclinations o! t he ascent trajoctoiv planes ate functions ol 
launch site latitude 0 t and azimuth I 7 only 

cos i ‘ cosO ; V sin (151 

The two darlv opportunities do drttei greatly . !n*we\er. in 
the light ascension ol the a$» ending mode 11 ot the orbit and 
in the length ot the travel sed in-plane ate, the ranee angle #1. 


sin ia ; il) 


sin \ sm 1‘j 


Ouadrant rules tor this equation involve lire observation that a 
negative cos places (it 7 - ii) into the second oi thud quad- 
laid, while the sign id sin (o 7 S2) determines the choice 

between them. 

I he lange angle ft is measuted in the menial ascent tiaiec- 
ton plane Itoni the lift -oil point at launch all the wav to the 
depaituic asvmptote duection. and can be compuied loi a 
given launch tune t t oi a,. - o 7 u 7 | and an azimuth 
ahead v known from l*q ( l M as follows 

cos - sin iS m \ sin y* + co*\ \ cosO y V cos (a, a ) 
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-90 -70 -50 -30 -10 10 30 50 70 90 

ASYMPTOTE DECLINATION , deg 

Fig. 13. Permissible regions of azimuth vs asymptote 
declination launch space for Cape Canaveral 


sin 0 


sin (or^ - or. ) X cos 8^ 


( 18 ) 


extent of range angle 0 can be anywhere between 
0 deg and 360 deg, so both cos 0 and sin 0 may be desired in 
its determination. The range angle 0 is related to the equa- 
torial plane angle, Aar = a m - cl l , discussed before. Even 
though the two angles are measured in different planes, they 
both represent the angular distance between launch and depar- 
ture asymptote, and hence they traverse the same number of 
quadrants. 

Figure 1 5 represents a generally applicable plot of central 
range angle 0 vs the departure asymptote declination and 
launch azimuth, computed using Eqs. (17) and (18) and a 
launch site latitude </> L = 28.3 (Cape Canaveral). The twin 
daily launch opportunities are again evident, showing the 
significant difference in available range angle when following a 
vertical, constant 8 M line. 

It is sometimes convenient to reverse the computational 
procedure and determine launch azimuth from known range 
angle 0 and t RLT > i.e., Aa = - a^, as follows: 



RIGHT ASCENSION a, deg 
(EARTH EQUATORIAL CELESTIAL LONGITUDE) 

Fig. 14. Typical launch geometry example in celestial (inertial) Mercator coordinates 
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DEPARTURE ASYMPTOTE DECLINATION 6,,, , deg 


* 

Fig. 15. Central range angle 0 between launch site and outgoing asymptote direction vs its declination and launch azimuth. 
Pair of typical example launch windows for 6^ * - 10 deg are shown by bold line segments 

(reproduced from Ref. 7). 


snl L 


cos 6. X sin (a. - <* t ) 
sin 0 


sin - cos 0 X sin 
C0S ‘‘7- sin 6 X cos <(> L 


with the lighting conditions at lift-off and consequently allows 

( 19) a lighting profile analysis along the entire ascent arc. 

The angle ZALS, displayed hi Fig. 16, is defined as the 
angle between the departure V„ vector and the Sun-to-Earth 

(20) direction vector It allows some judgment on available ascent 
lighting. 


Figure 16 displays a 3-dimensional spatial view of the same 
typical launch geometry example shown previously in map 
format in Fig. 14. The difference in available range angles as 
well as orientation of the trajectory planes for the two daily 
launch opportunities clearly stands out. In addition, the figure 
illustrates the relationship between the “first" and “second 
daily" launch windows, defined in asymptote-relative time. 
t KlT . as contrasted with “morning” or “night" Lunches, 
defined in launch-site-local solar time. The latter is associated 


The length of the range angle required exhibits a complex 
behavior- the first launch window of the example in Figs. 14 
and 16 offers a longer range angle than the second, but the 
second launch window opens up with a range angle so short 
that direct ascent into orbit is barely possible. Further launch 
delay shortens the range even further, forcing the acceptance 
of a very long coast (one full additional revolution in parking 
orbit) before transplanetary departure injection. A detailed 
analysis of required arc lengths for the various sub-arcs of the 
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Fig. 18. Angle from perigee to departure asymptote 


= is the true anomaly of the injection point, usually 
is near 0 deg, and can be computed by iteration 
using: 


tan?/ 


e u smv i 
1 + e H cos 


(23) 


7, = injection flight path angle above local horizontal, 
deg 

v n - eccentricity of the departure hyperbola: 


‘ ' 1 + 


C ’< X 'r 


(24) 


To make fcq. (22) balance, the parking orbit coast arc, 
6 coast , must P ic ^ U P any slack remaining, as shown in 
Fig. 17. A negative 0 coaST implies that the 1, -solution was 
too short-ranged. A direct ascent with positive injection true 
anomaly v t (i.e., upward climbing (light path angle, 7 ; , at 
injection) with attendant sizable gravity losses, may be accept- 
able. or even desirable (within limits) for such missions. 
Alternately, the other solution for I, . exhibiting the longer 
range angle 0. and thus a longer parking orbit coast. S COASr . 
should he implemented An extra i evolution in parking orbit 
may be a viable alternative. Other considerations, such as 
desire for a lightside launch and/oi injection, tracking ship 
location and booster impact constraints, may all play a signif- 


icant role in the ascent orbit selection. A limit on maximum 
coast duration allowed (fuel noil-off, battery life, guidance 
gyro drift, etc.) may also influence the long/short parking orbit 
decision. In principle, any number of additional parking orbit 
revolutions is permissible. Shuttle launches of interplanetary 
missions (e.g., Galileo) are in fact required to use such addi- 
tional orbits for cargo bay door opening and payload deploy- 
ment sequences. In such cases, however, the precessional 
effects of Earth’s oblateness upon the parking orbit, primarily 
the regression of the orbital plane, must be considered. 

4. Parking orbit regression. The average regression of the 
nodes (i.e.. the points of spacecraft passage through the 
equator plane) of a typical direct (prograde) circular parking 
orbit of 28.3-deg inclination with the Earth's equator, due to 
Earth's oblateness, amounts to about 0.46 deg of westward 
nodal motion per revolution and can be approximately com- 
puted from 


540° X r l X 7, X cos i 

S 2 = : , deg/revolutu»n (25) 

r 2 

(> 


where 

r = Earth equatorial surface radius. 6378 km 

- circular orbit radius, typically 6748 km for an 
orbital altitude of 370 km ( 200 nmi) 


16 




* ■ I 1 : ; 

, r n.- 


J 2 - 0.00108263 for Ear iii 

i = parking orbit inclinatioii. dc s * * computed for 
a given launch geometry fr* s 5 

- V (26) 

This, cmrection. multiple / .he orbital stay time of A' 
revolutions, must be c-onskl. »» 2 n; determining a biased launch 
time and, hence, the ngh» ascension of the launch site at 
lift-off: 

a, - u. + il X A' , deg (27) 

eff L 

S. Dogleg ascent. Planar ascent has been considered exclu- 
sively, thus far. Reasons for performing a gradual powered 
plane change maneuver during ascent may be many. Inability 
to launch in a required azimuth direction because of launch 
site constraints is the prime reason for desiring a dogleg ascent 
profile. Other reasons may have to do with bum strategies 
or intercept of an existing orbiter by the ascending spacecraft, 
especially if its inclination is less than the latitude of the 
launch site. Doglegs are usually accomplished by a sequence of 
out-of-plane yaw turns during first- and second-stage bum, 
optimized to minimize performance loss and cominenc* ig as 
soon as possible after the early . low -altitude, high aerodynamic 
pressure phase of flight is completed, or at ter the necessary 
lateral range angle offset has been achieved. 

By contrast, powered plane change maneuveis out of 
parking orbit or during transplanetary injection are much less 
efficient, as a much higher velocity vector must now be rotated 
through the same angle, but they may on occasion be opera- 
tionally preferable. 

As already discussed, a special geometric situation develops 
whenever the departure asymptote decimation magnitude 
exceeds the latitude ot the launch site, causing a "split azi- 
muth" daily launch window Figure 13 shows the effects of 
asymptote decimation and range satety constraints upon the 
launch problem As the absolute value ot declination increases, 
it eventually reaches the safety constraint on azimuth, prevent- 
ing any further planar launches The situation occurs mostly 
early and or late in the mission's departure launch period, and 
is frequently associated with dual launches, when month-long 
departure periods are desired 

The Shuttle-era Space Transportation System (STS), includ- 
ing contemplated upper stages, is capable of executing dogleg 
operations as well. These would, however, effectively reduce 
the launch vehicle's payload (or C\l capability, as they did on 
expendable launch vehicles ot the past. 


6. Tracking and orientation. As the spacecraft moves away 
front the Earth along the asymptote, it is seen at a nearly 
constant declination that of the departure asymptote. 
(DLA in the plotted data). The valrc ot D! A greatly affects 
tracking coverage by stations located at various latitudes; 
highest daily spacecraft elevations, and thus best reception, are 
enjoyed by stations whose latitude is closest to DLA. Orbit 
determination, using radio doppler data, is adversely affected 
by DLA's near zero degrees. 

The spacecraft orientation in the first few weeks is often 
determined by a compromise between communication i antenna 
pointing) and solar heating constraints. The Sun-spacecraft- 
Earth (SPE) angle, defined as the angle between the outgoing 
K-infinity vector and the Sun-to-Earth direction, is very 
useful and is presented in the plots under the acronym ZALS. 
It was defined in the discussion of Fig. 16. This quantity has 
many uses: 

(1) If the spacecraft is Sun-oriented, ZALS equals the 
Earth cone angle (CA), depicted in Fig. 19 (the cone 
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an jnulc between the vehicle's longitudinal -/axis ..mi 
the object dnectiuii). 

(2) The Sun-plase angle. <l\. (phase angle is the Sun-object - 
spacecraft angle) describes the state of the object's 
disk lighting: a fully lit disk is at zero phase. For the 
Earth (and Moon), several days after launch, the su.i- 
phase angle is: 

<h v - I M) - / A I s [IS) 

(3) The contour labeled /ALS ” W separates two cate- 
gories of tianster traiec tones those cail> departures 
that first cut inside hjrdfs orbit. thus stalling out at 
negative hehoccntuc true anoinalie. tor /ALS > K K) 
and those later one* that stait at positive tiue anoma- 
lies. heading out toward Jupitci and never expenemmg 
the incieased soldi heating at distances ot less tlian 
1 Al tor /ALS v oo* 

7. Post-launch spacecraft state. Alter the spacer rail has 
departed ft out the immediate vicinity ot Faith ti.c.leU the 
I at tli’s *phcic ot influence o| about 1-2 million km), it moves 
on a heliocentric conic whose initial conditions may be 
approximated as 


^ N< **/ 1 A/// * Vs*i ‘ ^ 

where K and V ot I jith aic evaluated Iron 1 an ephemeitN at 
lime ot miection and A/ lepiesents tune elapsed Miue then <m 
seconds) 1 he vector m I Ml 50 ^aitesian * ooidm.de van 
be constructed using \ (\t I \ I>1 A - \ and Kl A o 
m thsec components as tolious 

|l^ ■ COS ’ cOSiS^. ' 

r„ - sm \ \\\ i 

H. Orbital launch problem Orbital launch tiom the Shuttle, 
ttoni oilier elements ot the SIS -u trom anv tempmuiv or 
permanent oibital space station complexes, introduces eiituelv 
new concepts uito the l.ailh-depaituic problem Some ot 
the new -mstiamts. ah cad v mentioned, limit out ahluv 
to launch a given uileip. mctaiy mission llvsiowlv tegtessing 
space station oibit (see 1\] 25) genet all> does not contain 
the 1 intimtv vedoi reijuited at depattuie Dibit htetime 
oi other consultations ma> dictate a .pave station's orbital 
altitude that mav he too high tor an elticieui injection hum 
Innovative departure strategies jre beemmne to erne tee. 
attempting n- alleviate these problems a re*ent Sdeiue 


Applications. Inc. (SAM studj (Re!. u ) rnjiifj fi% cimi.j nt’ 
the techniques available, such as passive wait for natural 
alignment of die continuously regression .space station orbit 
plane (driven hy F.arth's oblatencss) witn the required 
' -infinity vector, or the utilization of 2- and 3- impulse man- 
euvers. seeking to perform spacecraft plane changes near the 
apogee of a phasing orbit where velocity is lowest and thus 
turning the orbit is easiest. These tvv-> approaches can he 
combined with each other, as well as wi*h other suitable 
maneuvers, such as. 

( 1 ) Deep space ptopulsivc burns fot orbit shaping and 
phasing. 

t2) (iravny assist tlyby. including Farth return AYFfiA, 

1 31 Aeiody namic turns at grazing perigees oi at inter- 
mediate planetary s»« mgbvs.aml 

(4) Multiple rev ol utit >n tniecMon burns, requiring several 
low . grazing passes, combined with apogee plane change 
maneuvers, etc. 

All u| these devices can be optimized to permit satisfactory 
orbital launches, as well as to achieve the most desirable condi- 
tions ji the final arrival body - In general, space launch advan- 
tages. such as on-oi bit assembly and checkout of pay loads and 
clustered multiple propulsion stages, or orbital construction of 
bulky and P agile subsystems (solar panels, sails, unieiiuas. 
radiator n. booms, c.c.i will, it is hoped, greatly outweigh the 
significant deep-space mission pe, .allies incurred because of 
the space station's inherent oibital orientation incompatibility 
with dcpaiturc requirement 

D. Planetary Arrival Synthesis 

I lie planctaiv arrival traject >iv design problem involves 
satisfying the project's engineering and science objectives at 
lire target body by duping the aitival trajectory m a citable 
maiinei. As these objectives mav be quite diverse, only io*:t 
illustrative scctianos shall be discussed i.i this section tlv'v. 
or bile! atmosphere probe, and to a small extent .ai.dei 
missions 

I Flyby trajectory design. In this mission mode, he arrival 
trajectory is not modified m ato deteimmistic wjv at the 
planet the oiigm.il aim point and amval Mine are chosen to 
sutistv trie largest number ot potential objectives, long he tote- 
band 

Ibis pto„ess involves the choice t-l jinval dati to ensure 
desirable characteristics. siuh as t l »e values ot the variables 
Mil*. b\P /APS. cU . presented in plotted toim in the data 
se,. ti«*n ot this \ olume 
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l>;\P. the planet (*«|ii.iI*mi.i) iUvhn.it ion. . ot the incoming 

as\ mptole. i c., ot tho 1 uitinin vectoi, psovuies the measiue 
ot the nuuunum possible inclination ot ll\h\ Its negative is 
also Known as the latitude ot leitical uup.u t (l VIK 

The magnitude ot I iiitimh , VIII 1 IV !, enables one to 
oontiol tiu* tl>h> fm ii .mi'll' Si between llu* meouuiig atul 
outgoing V bi vivlois In a suitable choice ot closest appioach 
t< ' U i.ulniN, / /( (mv l u* HU 

Av ISO . deg t 

wlieic /». the as\ mptole lull angle, is I omul liom 


1 1 

cos ft ~ (.«.U 

* I : / 


VHP also enables the designei to evaluate planet oconluc 
velocity. I*, at anv distance . /. on the llyh\ liypeibola 

/ ' /i t 

I l J k l \ . Kin x l.MI 

In the o.hcvc equations. ji ;i (oi (».W >. in the giavitalional paiam- 
etc i ot the auival bod\ . 
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Another pair of significant variables on which to base 
arrival date selection are ZAPS and ZAPE-the angles between 
V infinity and the phnet-to-Sun and -Far»b vectors respec- 
tively. These two angles represent the cone angle (CA) of the 
planet during the far-encounter phase for a Sun- or an Earth- 
oriented spacecraft, in that order. ZAPS also determines the 
phase angle. 4> s . of the planet's solar illumination, as seen by 
the spacecraft on its far-encounter approach leg to the planet: 

<f> 5 - 180 -ZAPS (35) 

Both the cone angle and the phase angle have already been 
defined and discussed in the Earth departure section above. 

The flyby itself is specified by the aim point chosen upon 
the arrival planet target plane. This plane, often referred to 
as the 5-plane, is a highly useful aim point design tool. It is a 
plane passed through the center of a celestial body normal to 
V^. the relative spacecraft incoming velocity vector at infinity. 
The incoming asymptote, i.e.. the straight line, zero-gravity 
extension of the V . -vector, penetrates the 5-planc at the aim 

point. This point, do d by the target vector B in the 5-plane, 

is often described by its two components B • T and B • R, 
where the axes T and R form an orthogonal set with V^. The 
T-axis is chosen to be parallel to a fundamental plane, usually 
the ecliptic (Fig. 2i) or alternatively, the planet's equator. The 
magnitude of B equals the semi-minor axis of the flyb\ hyper- 
bola. b. and can be related to the closest approach distance, 
also referred to as the periapse radius, r , by 

-2 2 \ ,/2 

"" -’) ■ km '• ,6> 
or 

kn ’ W7) 

The direction angle 0 of the 5-vector. B. measured in the 
target plane clockwise from the T-axis to tiie 5-vector position 
can easily be related to the inclination, i. of the flyby trajec- 
tory. provided that both (DAP) and the T-axis, from which 
0 is measured clockwise, are defined with respect to the same 
fundamental plane to which # he inclination is desired. For a 
system based on the planet equator (Fig. 22): 

cos i pFQ = cos 0 phv X cos b m ^ (38) 

which assumes that d PK y is computed with the T-axis parallel 
to the planet equator (i.e.. = V_ X POLEpj.yl, not the 

ecliptic, as is trequcntK assumed (T F(V = X POLE^.j ). 
Care must be taken to use the 0-auglc as defined and intended. 



B MISS PARAMETER. BJS 
(TARGET VECTOR) 

0 AIM POINT ORIENTATION 
S PARALLEL TO INCOMING ASYMPTOTE. V^, 

T PARALLEL TO ECLIPTIC PLANE 
ANDITOS 

R = Sx T 

Fig. 21. Definition of target or arrival 5-plane coordinates 

The two systems of 5-plane T-axis definition can be recon- 
ciled by a planar rotation. -A0. between the ecliptic T- and 
R-axes and the T- and R-axis orientations of the equator 
based system 

sin (a. - a p:p ) 

tan A0 = 1 — . £ ‘ v ; : 

cos X tan b pp - sin 6^ X cos (a w - o^) 

(39) 

where 

o KP and b Fp are right ascension and declination of the 
ecliptic pole in planet equatorial coordinates. For Jupiter, 
using constants in Section V : 

<x FP * 290.598 b Fp = 87.789, deg 

and b„ are RAP and DAP, the directions of incoming 
also in planet equatorial coordinates. 
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ECLIPTIC POLE (EP) 



« PLANE (I TO V.J 


J m " v « x p m. 

T /v c> ' v ~ x P /V-'U 

Ao ' 90 |(o A? - 90) - (o^ - 180)) a., - 

AO ARCTAN (SIN Aa/(COS 5^ x TAN S f / , SIN « x , x COS Aft)| 
Fig. 22. Two T axis dafinitions in the arrival plane 


The correction SO is applied to a 0 angle computed m the 
ecliptic system as follows (Fig. 22): 

°n V ",<•/ HOI 

I he ecliptic T. Rf.* i7 axes, how eve i. have to be totaled b> 
-SO (clockwise duectum is positive in the tfplane) to obtain 
planet eijiia tonal T. Rp^ cooidinate axes The //•magnitude 
of an aim point in eithei s> stem is the same 

Ihe pioiections of the Sundo-planet and Faiih-io-planet 
vectors into the li plane represent aim point loci of diametiic 
Sun and Faith occupations, lespectivelv, as defined in Fig. 23. 
The H plane 0 angles (\\ nh tespect to T fVi axis) of these van 
ables ate presented and labeled FTSP and F1TP, respective!) , 
in the plotted mission data. In addition to helping design oi 
else avoid diametiic occupations, these i|iiannties allow 
computation of phase angles. *b f , of the planet at the spaceciaft 


peiiapse. at the entix point ot a piobe, oi genet ally at an> 
position r (subsetipt A’ Sun. could be unplaced by F- Farth. 
if desired), (see Fig. 24): 

cos «l* v ■- -Ci's ^ \ ci's /\\' s - sin \ sin /Al\. 

\ l s r -0 sc ) ( 41 ) 

wheie 

r ,hc Jlin P oml angle n the H plane, must he w ith 
lespect to the same T . as F.TSP 

= amval mnge ingle fiom intimtv (a position tai 
out on the incoming asxmptote) to the point ot 
mteiest r (Fig 25) 

Ihe computatii'u of the anival lange angle, to the 
position ot the desited event depends on its t\pe. as follows 
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I . At flyby periapse. (0*, > ^ 0 deg): 


cos 0^ = cos(-pJ 


-1 



( 42 ) 


(r p is periapse radius). 

2. At given radius r . anywhere on the flyby trajectory 
(see Fig. 25): 


+ v r 03 ) 

f w should be computed from periapse equation. 
Eq. (42) v f at r can be obtained from < v f + 
v^, v f has negative values on the incoming branch): 



JS. At the entry point hiving a specified flight path angle 
V U ; ig- 25). 


Fig. 23. Definition of apuroach orientational coordinates 
ZAPS and ETSP. ZAPF. and ETEP 




+ r,. 


(45) 


«-PLANE 


ECLIPTIC 



TO SUN 

TRAJECTORY ARC 


tX'l. 


TERMINATOR 


r. RADIUS TO 
A POINT ON 
TRAJECTORY 
PHASE ANGLE. «l> PROJECTED ONTO 


UNIT SPHERE 
A ETSP - 180° - ft 


Fig. 24. Phase angle geometry at arrival planet 
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where v F is the true anomaly at entry, should always 
he negative, and can be computed if entry radius and 
altitude, r E = r SURF + h E and y t:% the entry angle, are 
known: 



whereas the fictitious periapse radius r p for the entry, 
to satisfy y h: at r t> is equal to 



(47) 


The B value corresponding to this entry point can be com- 
puted from Eq. (36), while 0 ,. de in the fFplane would 
depend on the desired ent*., latitude inclination. Eq. (38). or 
phase angle, i.q. (41). 

The general flyhy problem poses the least stringent con- 
straints on a planetary encounter mission, thus allowing 
optimization choices from a brge list *.f secondary parameters, 
such as satellite viewing and occultation. planetary fields and 
particle in situ measurements, special phase-angle effects, etc. 

A review of the plotted handbook variables, required in the 
phase-angle equation (Eq. 41). shows that the greatest magni- 
tude variations are experienced by the ZAPS angle, which is 
strongly flight-time dependent: the longer the trip, the smaller 
ZAPS. For low equatorial inclination, direct flyby orbits, this 
implies a stead\ move of the periapse towards the lit side and. 
eventually, to nearly subsolar periapsesfor long missions. This 
also implies that on such (lights the approach legs of the tra- 
jectory are facing the morning terminator or even the dark 
side, as trip time becomes longer, exhibiting large phase angles 
(recall that phase is the supplement of the ZAPS angle on 
the approach leg). This important variation is caused by a grad- 
ual shift of the incoming approach direction, as flight time in- 
creases. from the subsolar part of the Jovian leading henmpheie 
(in the sense of its orbital motion) to its antisolar part. 

The arrival time choice on a ver> fine scale may greatly 
depend on the desire to observe specific atmospheric, 'surface 
feat 111 es (e.g., the Great Red Spot) or to achieve close encoun- 
ters with specific satellites of the arrival planet. Passages 
through special satellite event zones, e.g.. flux lubes, wakes, 
get ven trie and 01 heliocentric occupations, require close 


control of arrival time. The number of satellites passed at 
various distances also depends on the time of planet C//4 . 
These fine adjustments do, however, demand arrival time 
accuracies substantially in excess of those provided by the 
computational algorithm used in this effort, which generated 
the subject data (accuracies of 1-5 min for events or 1-2 h 
for encounters would be reo 1 ired vs uncertainties of up to 
1 .5 days actually obtained with the rectilinear impact pseudo- 
state theorem). Numerically searched-in integrated trajec- 
tories. based on the information presented as a first guess 
input, are mandatory for such precision trajectory work. 

Preliminary design considerations for penetrating, grazing, 
or avoiding a host of planet-centered fields and particle struc- 
tures. such as magnetic fields, radiation belts, plasma tori, ling 
and debris structures, occupations by Sun. Earth, stars, or 
satellites, etc., can all be presented on specialized plots, e.g.. 
the A-plane. and do affect the choice of suitable aim point 
and arrival time. All of these studies require the propagation 
of a number of flyby trajectories. Adequate initial conditions 
for such efforts can be found in the handbook as: VHP (Kj 
and DAP (6^) already defined, as well as RAP ( 0 ^). the 
planet eqi atoria* right ascension of the incoming asymptote 
(i.e.. its east longitude from the ascending node of the planet's 
mean orbital plane on its mean equator, both of date). The 
designer's choice of the aim point vector, either as B and 0, or 
as cartesian B * T and B * R. complete the input set. Suitable 
programs generally exist to process this information. 

2. Capture orbit design. The capture problem usually 
involves the task of determining what kind of spacecraft orbit 
is most desired and the interconnected problem of how and at 
what cost su;h an orbit may be achieved. A scale of varying 
complexity may be associated with the effort envisioned an 
elliptical long period orbit with no specific orientation at the 
trivial end of the scale, through orbits of controlled or opti- 
mized lines of apsides (i.e.. periapse location), nodes, inclina- 
tion. or a safe perturbed orbital altitude. Satellite G/A-aided 
capture, followed by a satellite tour, involving multiple satel- 
lite G/A encounters on a number of revolutions, each designed 
to achieve specific goals, probably rates as the most complex 
capture orbit class Some orbits arc energetically very difficult 
to achic/e, such as close circular orbits, but all require signifi- 
cant expenditures of fuel. As maneuvers form the background 
to this subject a number of useful orbit design concepts shall 
be presented to enable even an unprepared user to experiment 
with the data presented. 

The simplest and most efficient mode of orbit infection is 
a coplanar hum at a common periapse of the arrival hyperbola 
and the resulting capture orbit (Fig. 26). The maneuver 41’ 
required is: 
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Si 





, km 


(4S) 


A plot of orbit mseition SI' u\] lined as a function ot /j, ami 
/* lusinji I tj. 50) is piesented in l : ig. ?**. Hie apoapse ladtusof 
such an oibit of given period would be 


The oibital peuod foi such an mbit, tequmnj: knowledge 
ot per lapse aiul apoapse tadii. r f% and r 4 . is 


/’ 



(4 l >) 


It' on flic olhei hand, a known oibit peuod /* (in seconds) 
isdesued. the expiession tot AT is 




An evaluation ot* l\|. (>0) (and l ; ig. - *) shows that lowest 
oibit mseition is obtained foi the lowest value i»i 'r the 
ionite »t peuod l\ and the lowest I of ai rival 

Ot some mteiest is inieetion into citcular captuic oibits. a 
special case of the coapsulal insertion pioblem. Ii can be 
shown (Ret. S) dial an optimal A I exists foi mseition into 
cap lute oibits of constant eccentiicitv. including «• 0. u\. 

cnculai oibits. which would lequiie a specific ladius. 



Fig. 26. Coapsidal and cotangantial capture orbit insertion geometries 
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Fig. 27. Coapsidal capture orbit insertion maneuver AV 
requirements for Jupiter (using Eq. 50) 


Au > p , one can solve for the hyperbolic periapse r CA and the 
burn radius r using selected values of true anomaly at hyper- 
bolic burn point v H and it« capture orbit equivalent 


•v = v „ + 


(SS) 


utilizing the following three equations (where hi and H stand 
for elliptic and hyperbolic, respectively): 


and 


sin v ff ca 






SMI I' / ,-2 v 

—v) 


CA 


km 


r = 




+ 1 


(56) 


(57 1 


r o> = 


2/i 

»» 

r 2 


(52) 


while the corresponding optimal value for A F would be 




(58) 


I 

uc 



(53) 


Frequently the orbital radius obtained by use of Eq, (52) is 
incompatible with practical injection aspects or with arrival 
planet science and engineeiing objectives. 


A more general coplanar mode of capture orbit insertion, 
requiring only tangent iality of the two trajectories at an 
arbitiaiy mainover point of udiusr common to both oibits. 
Fig 2h. teq uires a propulsive elToit of 


Al 



/ 2 W+ r r- 

7 r ^A +r r ) 


(54) 


It can be dearly seen that by performing the burn at periapse 
the substitution r = r brings us back to Eq. (48). 

The cotangential maneuver mode provides nonoptimal con- 
trol over the orientation of the major axis of the capture orbit. 
If it i\ desiied to rotate this line of apsides clockwise by 


The procedure of obtaining a solution to these equations is 
iterative. For a set of given values for r, , r . and an assumed 
Aco. a set of v H and i’^. the hyperbolic and elliptical bum 
point true anomalies which would satisfy Eqs. (55-58) can 
be found. This in turn leads to r, the maneuver point radial 
distance, and hence. AT (Eq. 54). A plot of AT cost for a set 
of consecutive A choices will provide the lowest AT value 
for this maneuver mode. 

For an optimal insertion into an oibit of an aibiti aiy niajoi 
axis orientation one must turn to the nioie genet al. still co- 
planar. but intersecting (i.e.. noutangeniial bum point) man- 
euvei (see Fig. 28). It provides sufficient flexibility to allow 
numerical optimization of Al’ with respect to apsidal rotation. 

Au> 

p 

A more appioptiale way to define apMdul orientation is ti> 
measuie the post -maneuver capture orbit periapse position 
angle with respect to a fixed direction, e.g.. a far encounter 
point on the incoming asymptote. -V^, thus defining a cap- 
ture oihit periapse range angle 
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Fig- 28 . Coapsidaf and intersecting capture orbit insertion gaometries 

where 


- Acj 


+ arc cos / \ 


(59) 


Taken Irom Ref. 8. tl>e expression for the intersecting burn 
•ills 


^ fUT 

'l' '.%J r’ J * Q 


Q = | J^ Xr r*'cA r i I + 

* >/£ ^ < r " 'or > + »’i I ' 2 - Ta )T^;, (T - ,)" J 

(60) 

r ~ planet-centered radius at burn, r ’* r .'r . 
r i a,,J r r ~ apoapse and penapsc radii ol capture ellipse 


n 
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r CM = closest approach radius of flyby hyperbola 

£ * flag: 6 = +1 if injertinn occurs on same leg 

(inbound or outbound) of both hyperbola 
and capture ellipse, 6 - -1 if not. 


It should be pointed out that Eq. (57) and (58) still apply 
in the intersecting insertion case, while Eq. (56) does not 
(as it assumes orbit tangency at burn point). 

The evaluation of intersecting orbit insertion is more 
straightforward than it was for the cotangential case. By assum- 
ing K, and orbit size (e.g., V m s 6 km/s, = 4 RJ, P = 200 
days, similar to a Galileo capture orbit) and stepping through a 
set of values for v £ . the capture orbit burn point true anomaly, 
one obtains, using Eqs. (57-60), a family of curves, one for 
each value of ^CA- the hyperbolic closest approach distance. 
As shown ip. Fig. 20, the envelope of these curves provides the 
optimal insertion burn AV for any value of apsidal rotation 
AU", desired. The plot also shows clearly that cotangential and 
apsidal insertion burns arc energetically inferior to burns on 
the envelope locus. For the same assumed capture orbit, a 
family of optimal insertion envelopes, for a range of values 
of arrival is presented in Fig. 30. The smallest value. 
F, = 4.458 km/s, represents the remaining pre-insertion 
energy after an lo encounter, typical of Galileo mission 
strategy. 


The location of periapse with respect to the subsolar point 
is of extreme importance to many mission objectives. It can be 
controlled by choice of departure and arrival dates, by AV 
expenditure at capture orbit insertion, by an aerodynamic 
maneuver during aerobraking, by depending on the planet's 
motion around the Sun to move the suhsolar point in a manner 
optimizing orbital science, or by using natural perturbations 
and making a judicious choice of orbit size, equatorial inclina- 
tion, i, and initial argument of periapsis. w 0 , such as to cause 
regression of the node. 12, and the advance of periapsis, cb. 
both due to oblateness to mo/e the orbit in a desired manner 
or at a specific rate. Maintenance of Sun-synchronism could 
provide constant lighting phase angle at periapse. etc., by some 
or all of these techniques. For an elliptical capture orbit 
(Fig. 31). 


ii 



cos / X 


180 . , 
, deg/s 


( 61 ) 


? Rl»J, 

to x -V X — : ~ (2 -(5/2) sin 2 i) 


X 


(62) 


n 



, mean orbital motion, rad/s 


(63) 


A P 

a = — - — , semi-major axis of ellipitical orbit, km (64) 
2 r r 

9 A 9 P 

p = — - — ,semi-latus rectum of elliptical orbit, km (65) 

r A + r /» 

Rg ~ Jupiter's equatorial surface radius, km 

J 2 = Jupiter's oblatencss coefficient (for values see 
Section V on constants). 

For the example orbit (4 X 272.4 RJ) used in Figs. 28 
and 29, if near equatorial, the regression cf the node and 
advance of periapse, computed using Eqs. (61 -65), would 
amount to a negligible -0.23 and +0.46 deg/year, respectively. 
For contrast, for a grazing (a = 71.398 km) low-inclination, 
near-circular orbit, the regression of the node would race along 
at -64.5 deg/day. while periapse would advance at 129 d?g/ 
day. i.e.. it would take less than 3 days for the line of apsides 
to do a complete turn. 

It should be noted that 12-0 occurs for i = 90 deg, while 
cb = 0 is found for i = 63.435 deg. Sun-synchronism of the 
node is achieved by retrograde polar orbits, e.g., a 5 RJ cir- 
cular orbit should be inclined 1 1 1 .08 deg. 

3. Entry probe and lander trajectory design. Entry tra- 
jectory design is on one hand concerned with maintenance of 
acceptable probe entry angles and low relative velocity with 
respect to the rotating atmosphere. On the other hand, the 
geometric relationship of entry point, subsolar point and 
Earth (or relay spacecraft) is of paramount importance. 

Lighting during entry and descent is often considered the 
primary problem to he resolved. As detailed in the flyby and 
orbital sections above, the choice of trip time affects the value 
of the ZAPS angle which in turn moves the entry point tor 
longer missions closer to the subsolar point and even beyond, 
towards the morning terminator. 

Landers or balloons, regardless of deceleration mode, 
prefer the morning terminator entry point which provides a 
better chance for vapor-humidity experiments, and allows a 
longer daylight interval for operations following arrival. 

The radio-link problem, allowing data flow directly to 
Earth, oi via another spacecraft in a relay role, is very complex. 
It could require studies of the Earth phase angle at the entry 
locations or alternately, it could require detailed parametric 


where 
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Fig. 30. Minimum AV required for insertion into Jupiier capture orbit of given apsidal orientation 
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Fig. 31 . General satellite orbit paramatara and processional motion 
dua to oblatanaaa coefficient J 2 (from Raf. 8) 

studies involving relative motions of probe and relax spacecutt 
throughout probe enir> and its following slow descent. 

Balloon missions could also involve consideration of a 
vanet\ of wind dntt models, and thus. ate even more complex 
a , tar as the communications problem with the Faith oi the 
spacecraft is concerned 

E. Launch Strategy Construction 

I he omimi amts and desires. bricflx discussed abo\e. mas he 
displavcd oil the mission spate launJi amval da> plot as 
being limited bv the vontom houndaiies of /.. the dates 
1)1 A. VHP. /AP eu thus displaxtng the allowable lauiuli 

spa^e 

Within this lauiuli spa*.e a piefened da> b> da\ lauiuli 
strategy must he specified in a^woidaiue with piev ailing 
obtest tv cs Ihe simplest lauiuli sliategx .otten used to maintain 
a constant amval date at the Saieel planet icsulls in datlv 
lauiuh punts on a hon/ontal lute tnun lettmost tonghtmosl 
maximum allowable ( J boundaiv toi that amval date 
Such a stiateev makes use ot the lav r that most amval cliaiac* 
teristics max siav neailv constant acn*ss the launch space 
lighting and satellite positions in this case are fixed, thus 
allow me a similar eruounter satellite ti \. oi satellite tour 


A different choice of strategy could be to follow a contour 

lit #%♦ Cittna oli-ir'iof onelio uioli >i% ITI \ .u / I 

VIHMMS «^M « M«V>« « M ll . ‘/MV W«M»ISt 


also follow the minimum value locus of a parameter, e.c., 
C,/. (i.e . the houndarv between C lass I and 2 within Type I 
or Iff toi each launch dale, throughout the launch space. 


1 uiida:;iciilall> different is a launch si i a teg) toi a dual oi 
multiple spacecraft mission, involving more than one lauiuli. 
either tit winch rna> possibly puisue divergent objectives As 
an example Fig 32 shows the Voyagei . 1 and 2 lauiuh stiat- 
egy, plotted on an Farth departure vs Saturn arrival date plot. 
A 14-day pad turnaround scparc n between , **unches was 
to be maintained, a 10-day opportunity was to be available 
lor each launch, and the two spacecraft lud substantial!) 
different oh|ectives at Jupiter and Saturn one was to he 
lo-mtensive and a close J.ipitei flyby, to be followed by a 
close Titan encounter at Saturn, and the other was (ianymede- 
and or Callisto-intensive, a dist •!;» Jupiter tlybv.as a safety 
precaution agatns: Joviar radiation damage, aimed to continue 
past Saturn to Uranus and Neptune. Her:, even the spacecraft 
departure order was reversed by the strateg) within the 
launch space. 


Launch strategies toi oibita! departures from a space station 
m a specific orbit promise to introduce new dimensions into 
mission planning and design. New concepts are beginning to 
emerge on this subject e.g.. Refs. 9 and 10. 


O • • 



Fig. 32. Voyager 1MJS77) trajectory space and launch strategy 
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IV. Description of Trajectory Characteristics 
Data 

A. General 

The data represent trajectory peiformance information 
plotted in the departure date vs arrival date space, thus defin- 
ing all possible direct ballistic transfer trajectories between the 
two bodies within the time span considered for each oppor- 
tunity. Twelve individual parameters are contour-plotted. The 
first, C 3 /,, is plotted bold on a Time of Flight (TFL) back- 
ground; the remaining ten variables are plotted with hold con- 
touring on a faint C 3 Z background. Fleven plots arc presented 
for each of twenty mission opportunities between 19*5 an j 
2005. 

The individual plots are labeled in the upper outei corner 
by bold logos displaying an acronym of the variable plotted, 
the mission's departure year, and a symbol of the target planet. 
These permit a quick and fail-safe location of desired informa- 
tion. 

B. Definition of Departure Varabtes 

CjA: Farth departure energy (km 2 /s 2 ); same as the 
square of dt^ .mure hyperbolic excess velocity 
l l = CyL = v) - 2 where 

l' f = conic injection velocity (km/s). 

R f = R s + h /t injection radius (km), sum of 
surface radius R Spi A v/ . j injection 
altitude h r where Ksf. A R nf refers to 
Earth’s surface radius. (For valu*, see 
Section V on constants.) 

^ 9 gravitational constant times mass of *he 
launch body (for values, refer to Section V 
on constants) 

CjL must be equal to or exceeded by the launch 
vehicle capabilities. 

DLA; . geocentric declination (vs mean Farth equator 
of 1950.0) of the departure vector. May im- 
pose launch constraints (deg). 

RLA: a geocentric right ascension (v\ mean Farth 
equator and equinox of l l >50.0) of the departure 
\\ vector. Can be used with C 3 /. and DLA to 
compute a heliocentric initial state for traiectory 
analysis (deg) 

ZALS: Angle between depature vector and Sun-Faith 

vector. Equivalent to Farth-prohc-Sun angle several 
days out (deg). 

C. Definition of Arrival Variables 

V 1 IP I - 4 . planetocentnc arrival hvneibolu excess 
veloc:t> or l -infinity (km s). the magnitude of the 


vector obtained by vectorial subtraction of the 
heliocentric planetary orbital velocity from the 
spacecraft arrival heliocentric velocity. It repre- 
sents planet-relative velocity a* -reat distance from 
target planet, at beginning o. for encounter. Can 
be used to compute spacecraft velocity at any 
point r of flyby, including CIA (periapse) dis- 
tance r p : 



where 

Vpfjrprrnk = Rotational parameter GM 
SYS it: \n of the arrival planet system 

Jupiter plus all satellites. (For 
values, refer to Section V on 
constants.) 

DAP JL 4 , planetocentric declination (vs mean plane* 
equator of date) of arrival l’ m vector. Defines 
lowest possible flyby/orbiter equatorial inclination 
(deg). 

RAP o^, planetocentric right ascension (vs mean 
planet equator and equinox of date, i.e., RAP is 
measured in the planet equator plane from ascend- 
ing node of the planet's mean orbit plane on the 
planetary equator, both of date) Can be used 
together with VHP and DAP to compute an 
initial flyby trajectory state, bu ? requires & planc 
aim point informarion, e.g.. B and 0 (deg). 

ZAPS Angle between arrival l „ vector and the arrival 
planet-to-Sun vector. Fquivalent to planet-probe- 
Sun angle at far encounter; for subsolar impact 
would be equal to |K0 deg. Can be used with 
FTSP, VHP. DAP, and 0 to determine solar phase 
angle al periapse. entry, etc. (deg) 

ZAPF Angle between arrival I’* vector and the planet-to- 
Farth vector. Fquivalent to planet-probe- Far rii 
angle at far encounter (deg). 

FTSP Angle in arrival 6-plane, measured trom T-axis*, 
clockwise to projection of Sun-to-planet vector. 
Equivalent to solar »>ccultation region centerline 
direction in 6-plane (deg). 

FTFP Angle in arrival 6-plane, measured from T-axis , 
clockwise. to projection of harth-to-planet vector. 
Equivalent to Farth occultation region centerline 
direction in 6 plane (deg). 

•fcIS? and hTKP plot* are b *icd on T defined * being parallel to 
ciltptK plane ik , 1 text tor expUrutkiiil 
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V. Table of Constants 

Coi.'liinls list'll itt genei.ite (lu* ttitnini.ilion piesenled air 
stiiuinan/rd in this section. 

A. Sun 

i;.m n: kin’s-’ 

Vki .. / k "' 

B. Larth/Moon Systam 

UM M,? > 1 kn, < s- ’ 



*-.v. 

Knv ,> ‘ 

Suil.ue 
K 4(1 ms. 
km 

Mean 

Oil'll 

tUlnis.**’ 

kin 

Period. 
d«a> .> 

Jupitei 

i:6 (iHfr 4b: 

71.WS 0 


0 41 V.1R 

Kiiiftis 

SO \4 

ik:o 

4:t67| 

i ’bQi jk: 

1 IIIOp.1 

MOb 

\s.w 

r.7tW88 

I5MISIQ 

t«an\ turtle 

^S8s 

.'MIR 

1070)38 

7 IS4SSJ0 

t.llllvto 

7 i v 

:44.S 

188 '5 71. 

lb bRQOIS 


(•H'omputcd liom period ami Jupiter (»’M. miiiuled 
(•*11 ot System 111 Jupitei 10 tat ion ute 


<;\t, tKnl ws 1,00.44X0. •* km’ s-’ 

7, 0.00 1 OX .'(O 

Vk, „ v i*.*' , X 140 km 

C. Jupiter Systam 

'I .'MX \M km’ s : 

- 0.014 T .U00X 

'.n miK 


Direction of the Jovian planetary equatorial north pole 
(in l ; arlh Mean l-'quatoi of 1 ‘>50.0 coordinates): 
n ; , - .'t.X.OOl'W. ileit 

f> r - (,4.5040‘>. den 

0. Sources 

The constants represent the 1)1-1 IX planetary ephemeris 
( 1^X1 ) and the Voya(ter-2 Jupiter encounter trajectory* 
reconstruction data ( 107*)). 
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■52600.5 
-52500.5 
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■52300.5 
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2000 


Opportunity 


ENERGY MINIMA 


DEPARTURE ARRIVAL 

VALUE TYPE (YEAR/MONTH/DAY) (YEAR/MONTH/DAY) 
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84.783 
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